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Abstract 

In the primordial universe, fields with mass much larger than the mass-scale of the event- 
horizon (such as the Hubble parameter in inflation) exist ubiquitously, and can be excited 
from time to time and oscillate quickly around their minima. These excitations can induce 
specific patterns in density perturbations, which record the time dependence of the scale 
factor of the primordial universe, thus provide direct evidence for the inflation paradigm 
or its alternatives. Such effects are conventionally averaged out in theoretical and data 
analyses, but can be accessible for experiments targeting on density perturbations with high 
multipoles. 



1 Introduction 



The inflation [IHS], as the leading candidate paradigm for the primordial universe, has re- 
ceived strong support from observations of cosmic microwave background (CMB) and large 
scale structure (LSS) [1] . The simplest inflationary scenario not only explains the homogene- 
ity and isotropy of the Universe, but also predicts that the density perturbations seeding 
the large scale structure are generated at superhorizon scales, and are approximately scale- 
invariant, Gaussian and adiabatic, all of which are verified by experiments to some extent. 

Nonetheless, based on current observations, ambiguities and degeneracies still exit in 
terms of model-building. The specific inflation models remain illusive; in addition, there 
may be alternatives to inflation that have the same consequences on observables that we 
have been able to measure so far. To make further progress, there are at least two important 
questions. The first is how to find evidence that would unambiguously distinguish inflation 
from its alternatives. The second question is, if inflation is the correct paradigm, how to 
predict and measure new observables that will pin down the microscopic details. Similar 
questions apply to the alternative paradigms. 

Because all viable models have to satisfy the current observations that the two-point 
correlation function (power spectrum) of the density perturbations is approximately scale- 
invariant, for the purpose of the second question, studying small deviations from the scale- 
invariance and measurable higher-point correlation functions (non-Gaussianities) become 
very important. For example, for inflation models, properties of primordial non-Gaussianities 
can be classified [5], and if measurable, provide evidence for the interaction terms in the 
Lagrangian. Similar classification may be established for each of the alternative paradigms. 

However the first question still remains unanswered from this line of research. For ex- 
ample, assuming single field models and imposing the condition that the power spectrum is 
approximately scale invariant, the consequences on non-Gaussianities can be systematically 
worked out along the line of [6l4T0] , for either inflation or alternatives. The inflation models 
still predict approximately scale-invariant non-Gaussianities, while the attractor alternative 
paradigms |TTHT5] predict non-scale-invariant ones. However, as soon as we step away from 
this subset and consider multifield models, this sharp distinction will be lost. For exam- 
ple, for inflation models, multiple fields introduce various isocurvature modes that can have 
scale-dependent couplings to the curvature mode. Non-Gaussianities can be easily made 
non-scale-invariant if they are transferred from the isocurvature modes, for instance, when 
inflaton makes a non-constant turn in its trajectory in models of the type P^[T7] . Reversely, 
scale-invariant non-Gaussianities may be achieved in multifield or non-attractor single field 
non- inflationary models [T8|[T9] . In short, given general inflation scenarios, there is no generic 
prediction on how these non-Gaussianities should depend on scales. 

So far the primordial tensor mode is regarded as the only possible solution regarding to 
the first question. The tensor modes from inflation models are approximately scale-invariant 
with a red tilt; for some models they are observable. Typical alternatives such as the cyclic 
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model PUIET] or string gas cosmology [22p^ predict either non-observable tensor fluctuations 
or observable ones with blue-tilt. However, there are some important caveats for the tensor 
modes to achieve the goal unambiguously. Firstly, if we consider more general alternatives, 
scale-invariant and observable tensor modes are possible. The equation of motion obeyed by 
each polarization component of the tensor modes is the same as that by the massless scalar. 
So the tensor modes can be scale-invariant even in non-inflationary spacetime, just as the 
scalar. For it to be observable, we only need a large Hubble parameter. Scenarios of matter 
contraction [2ll[25] with large Hubble parameter are such explicit examples. Secondly, even 
for inflation, tensor modes are not guaranteed to be observable. While the best sensitivity 
for the tensor-to-scalar ratio achievable by experiments in the near future is Ar ~ O{10~^), 
the inflation models predict anywhere between r ~ O{10~^), for large fleld models, and 
r ~ O{10~^^), for small fleld models with TeV-scale reheating energy. 

So it is very important to search for complimentary properties in the density pertur- 
bations that can serve as a model-independent general distinguisher between inflation and 
alternatives. This is the main purpose of this paper. 

Before proceed, we would like to make a comment on the types of models we investigate. 
Arguably, inflation remains as the best available paradigm for the primordial universe. Its 
generic predictions naturally flt the data and its microscopic origin in term of fundamental 
theory is promising. Nonetheless, such opinions may be subject to personal taste; they are 
model-dependent and may even evolve with time. A more uncontroversial standard will 
be in terms of experimental data, and to ask what we can learn given the data by reverse 
engineering. So in this paper we will not discuss the important UV completion and model 
building aspects of the non-inflationary backgrounds. For the same reason, we will also not 
discuss which alternative of inflation is more natural than the others, for example, between 
expansion and contraction, attractor and non-attractor scenario. 

2 Bunch-Davies vacuum and resonance mechanism 

In nearly all models of primordial universe, the quantum fluctuations start their life in a 
vacuum that is mostly Bunch-Davies (BD). These fluctuations later exit the event horizon 
and become the seeds for the large scale structure. This applies to both inflationary and 
non-inflationary scenario, expansion and contraction universe, attractor and non-attractor 
evolution, single fleld and multifleld model, curvaton and isocurvaton modes. 

For example, consider the fluctuations of an effectively massless scalar fleld, S(j){'x,t), in 
a general time- dependent background with scale factor a{t), 



The conformal time r is deflned as dr = dt/a, and we will use dot to denote the derivative 
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with respective to t and prime to r. The event horizoE0 in physical coordinates is therefore 
|ar|. A quantum fluctuation with comoving momentum k is within the event horizon if 
k > 1/\t\. In this limit, the equation of motion for the fluctuations approaches that in the 
Minkowski spacetime limit. Along with the quantization condition, 

a^5(j)S(j)* - c.c. = i , (2.2) 

the mode function in the subhorizon limit becomes 

6(P -^=e-'^^ . (2.3) 
aV2k 



We have chosen the positive-energy mode, which corresponds to the ground state of the 
Minkowski spacetime, to be the BD vacuum. The effect of the background time-dependence 
is incorporated adiabatically in (12. 3p . The most important and universal property of (12. 3p 
is the oscillatory factor e~'^^'^ . Various prefactors depend on whether the form of Lagrangian 
(12.11) is canonical. 

To give explicit examples of the time-dependent backgrounds, we take the scale factor to 
be of the general power-law, 

a{t) = a{to){t/t^Y . (2.4) 

Because we require that the quantum fluctuations exit the event horizon, for p > 1 we need 
an expansion phase, so t runs from to +oo; for < p < 1 we need a contraction phase, so 
t runs from — oo to 0; for p < 0, we again need an expansion phase, so t runs from — oo to 0. 
The conformal time r is related to t by ar = t/{l—p)^ and r always runs from — oo to 0. For 
example, p > 1 corresponds to the inflation [l-3j, p = 2/3 the matter contraction pH[25]. 
p = 1/3 the pre-big-bang [26l[27], < p ^ 1 the ekpyrotic (slowly contracting) phase [20] . 
and — 1 ^ p < the slowly expanding phase [28j. 

To directly probe the universal BD vacuum, we need a high energy probe with wavelength 
much shorter than the event horizon. This can be achieved by introducing a small but highly 
oscillatory component in the background evolution p9]. Such a component resonates with 
the vacuum component which has the same physical wavelength. Because the BD vacuum 
has time-dependence, different momentum modes get resonated at different time. This effect 
is formulated in terms of the following integral, 

dTB{t)e-'^^ + c.c. . (2.5) 
The factor e"'^"^ in the integrand is the universal BD oscillatory component, and K is some 



^Here the event horizon is defined to be the maximum distance at t by which two points are separated 
but can still communicate with each other from t to tend- So it is a{t) dt/a = — ar, where Tend is set to 
0. 
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comoving momentum. The factor B{t) denotes the high energy probe mode we introduce. 
The integrand resonates when the two factors have the same frequency. Different momentum 
modes resonate one by one with the background, and in the mean while the phase of the 
background repeats due to oscillation. If we regard the repeated oscillation in B{t) as a 
clock, the time- dependence of the scale factor is translated into the ^-dependence of the 
integral through the resonance mechanism. For example if we take B{t) to be a periodic 
clock with frequency cu, B{t) ~ e**^*, f l2.5p becomes proportional tc§| 



~ sm 



— + phase 



p — 1 Hq \kr 



(2.6) 



where the "phase" denotes a i^- independent constant, kr is the mode that resonates at to; 
and Hq is the Hubble parameter at to- As we can see, the time-dependence of the scale 
factor is encoded inside the square bracket of (12.61) as a function of i^-modes. For power-law 
background and periodic resonance, this function is the inverse function of the scale factor. 
If we take the exponential inflation limit p ^ 1 and study a range of modes AK satisfying 
\n{AK/kr) < p, we have {K/Kf/P 1 + (1/p) \Yi{K/kr). So ([M]) goes to 



sm 



-— m - — h phase 
H kr 



(2.7) 



This is the leading resonance form found by Chen, Easther and Lim (CEL) for inflation [29] . 
As we can see, the CEL form is a special limit of the general resonance forms. In retrospect, 
the reason the argument of the sinusoidal function is proportional to In K is that this is the 
inverse function of the exponential function in inflationary scale factor. 

The distinctive oscillatory running behavior in the above resonant forms will not be 
changed by curvaton-isocurvaton couplings in multifield evolution. Any effect that also 
oscillates faster than the horizon time-scale generates additional resonance forms that su- 
perimpose onto each other. Any effect that varies much slower can only change the overall 
envelop of the resonance forms, by either changing their overall sizes or introducing scale 
dependent modulations. This latter modulation can also be informative as we will see in 
more details later. But similar to the scale-dependence of non-oscillatory correlation func- 
tions that we mentioned in Introduction, these scale-dependence can be rather arbitrary in 
multifield models, so much less robust than the resonant running. 

However, in terms of reverse engineering, we should also consider the possibility of non- 
periodic background oscillation components. A non-periodic background oscillation may 
cause resonance in a non-inflationary background, and conspire to give the same CEL form. 
For example, for arbitrary power law behavior (12.41) . we may engineer a background oscilla- 
tion component to be of the form B{t) ~ e^aMt/ta) gQ \]iqX the resulting resonance behavior 



^Using J^^dx e^^'^^^ ~ —e^^^^^'^y^e^^" / ^/±f^ , for positive/negative where the subscript de- 
notes the resonant point /'(a;*) = 0. Away from this point, the integrand e'-^^^^ is oscillating rapidly. 
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Ill - — h phase 

— 1 kr 
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which is the same as i\2.7\i . Therefore it becomes very important to search for standard 
clocks in physical systems. Such a clock should generate repeated perturbations with known 
time dependence, although not necessarily periodic. They should also be associated with a 
set of specific patterns that can be identified in observations. 



3 Spectator massive fields as standard clock 

Massive fields with mass much larger than the horizon mass-scale 1/ 1 ar | exist ubiquitously in 
models of primordial universejf] Even when we think of single field models, in a UV completed 
context, what we have in mind is really models with many massive modes. The single field 
model is obtained as the low energy limit where the energy scale is comparable to or smaller 
than l/|ar|, after these massive modes are integrated out. This is a good approximation 
even if the massive modes get excited classically and oscillate around its minimum. But 
for our purpose, these oscillations are a good candidate for the physical clock that we are 
looking for. So let us look at more details of the classical behavior of a massive particle a in 
the power-law background. 
The equation of motion is 



a + 3H& + mta = 



(3.1) 



where the Hubble parameter H = p/t. The solution is given in terms of Bessel functions. 
The asymptotic behavior of these Bessel functions at the limit rricrt ^ p^ is given in terms 
of sinusoidal functions, and we use these to approximate the oscillatory behavior of a, 



a 



-3p/2 



smim^t + a) + 



cos(mo-t + a) 



(3.2) 



where a is a phase, and a a is the initial oscillation amplitude at t = to- 
induce an oscillatory component to the Hubble parameter H, because 



Such oscillations 



1 1 

2 2 



other fields . 



(3.3) 



The leading term on the right hand side of (13. 3p does not oscillate in time because the 
energy is converting between kinetic and potential energy back and forth and conserved in 
the leading order. The oscillatory component for H, which we denote as -ffosd, comes from 



^For \p\ > 1, l/|ar| « \H\; for \p\ < 1, l/\aT\ > \H\ 
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Figure 1: A turning trajectory that excites the oscillation of massive fields. Dashed line 
indicates the potential valley. The massive field tends to settle down in the valley along the 
incoming and outgoing straight lines. But during the turning, the centrifugal force makes it 
deviate from the minimum. This induces the small oscillation. 

the subleading terms. Using (13. 2p . we get 



The next question is how these massive fields can get excited. There are many possibili- 
ties. As we have mentioned, even for single field models, we imagine a multifield configuration 
in which the effective single field trajectory turns from time to time depending on how the 
massive directions are lifted. During turning, the light mode and massive mode couple, so 
part of the energy can be released to excite the massive mode (Fig. [1]). The resulting oscil- 
lation ( 13. 2 p has a very high frequency. It can be averaged out in most cases, but not for our 
purpose. As we will see later in a more explicit example, even a tiny fraction of the energy 
transferred in this process can excite a large observable effect. More generally, massive fields 
may be excited classically by any sharp physical process, including the turning trajectory, 
sharp feature, particle creation and etc. 




(3.4) 



This in turn induces the oscillatory components for the parameters e = —H/H"^ and rj 
e/{He). Again we use the subscript "osci" to denote their oscillatory components. 




(3.5) 



(3.6) 
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4 Model and formalism 



We shall investigate two closely related processes and their observational signatures sepa- 
rately. One is the sharp feature that excites the massive fields. Another is the resonance 
phenomena induced by the oscillation of the excited massive fields. 

We use a two-field model in the general power-law background as an example. The sharp 
feature happens at to. Before to, the massive particle stays at minima and we consider the 
single field model 



-9 



(4.1) 



After to, we consider the two-field model 

1 



Co 



-9 



(4.2) 



Around to^ the a field is excited by some sharp process. For example, the field makes a 
turn (Fig. [T]). Note that after the turn, the two fields are still decoupled if it were not for the 
gravity. More complicated couplings are of course possible, but this minimum case is most 
general. Because the a- field is now considered as a spectator (except around to), to study 
the perturbation theory, it is important that we choose the following uniform-c/) gauge [17] . 
in which the scalar perturbation ^ corresponds to the conserved scalar degree of freedom in 
single field model. 



-N^df + h,j{dx' + Wdt){dx^ + N^dt) 



h 



1-3 



a e 



<Jt1 5 



, CT = cro(^) + ^cr(x, t) 



(4.3) 
(4.4) 



We use Maldacena's method [B] of the ADM formalism to expand the action. In this paper, 
we will only be interested in the correlation functions of so we ignore the perturbation Sa. 
The effect of a comes in because its zero- mode evolution o"o(t) perturbs the time-dependent 
couplings in the perturbative expansion. 

We separate the Hamiltonian in the perturbation theory as follows. 



aV-«eo(aC)' , 

-a^AeC' + aAe{dCf + C(Ae' 



(4.5) 
(4.6) 

(4.7) 



We have also separated e into the unperturbed part and the perturbed part due to features, 
e = eo + Ae, and in "Hg listed the only term important for this paper. The reason that this 
term is important is similar to that given in [29l[30] for infiation. Namely, the coupling in 
this term contains the highest time-derivative and becomes large in presence of features. 



7 



(b) 



Figure 2: Examples of Feynman diagrams used to perturbatively compute the power spectra 
and bispectra in feature models. 

Treating "Hg ^3 the interaction Hamiltonian, we can perturbatively compute the 
power spectrum and bispectrum using the in-in formalism, 



iCit)) ^ (o| 



T exp i / adrd X "H 



at) 



T exp 



adrd X W 



|0), (4.8) 



where the integration of r runs from —00 to 0. See |S] for a review of this formalism and 
methods for such computations. For example, the leading correction to the power spectrum 
is given by the diagram Fig. [2]^a), where the two-point vertex corresponds to (14.61) . 



1-0 

A(CkiCk.) = (OK / adrd^'x [^2^ Ck,Ck.]|0) . 



(4.9) 



Using the definition for the power spectrum P^, 



we get 



P 



CO 



(Ck,Ck.) = ^(27r)^5^(ki + k2) 



2i I dr a^Ae (m'^^ — ^i^^J + c.c. 
J —00 



(4.10) 



(4.11) 



where P^o is the power spectrum in the absence of features, and Uk{t) is the Fourier transform 
of C(x, t). The leading bispectrum is given by the diagram Fig. [2]^b), in which the three-point 
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vertex corresponds to (14.71) . 



•0 



(Ck.Ck.Cka) = {0\^ 



adrd^ 



(4.12) 



— CO 



(27r)^(5^(^ ki) + 2 perm. + c.c. . 




dul 



dr 



X 



(4.13) 



For simplicity, we will quote the bispectrum in terms of S{ki, k2, k^) according to the following 
definition [5], 



Using this method, other types of correlation functions in feature models can also be 
computed systematically. For example, the leading effect from the non-BD correction on the 
bispectrum [3l] corresponds to the diagram Fig. [2|^c). In this paper, we will only compute 
the two diagrams (14. lip and (I4.13p . which are the leading terms in this model. 

In inflationary and especially non-infiationary scenarios, there are variety of ways, in- 
volving one or more fields, to produce the leading scale-invariant power spectrum P^o iii 
absence of features. In this paper we do not concern how this is produced. We are interested 
in computing the resonance effects induced by massive fields, as corrections to the power 
spectrum and as the leading bispectrum. 

The above formalism applies to cases with arbitrary scale factor a{t). For the special 
case of inflation, different kinds of feature models have been studied in the past, including 
various effects from sharp features [32H3HI301I351I36] . periodic features [29l|371|3ll[3H], and 
massive particles |39tilT] . The main point of this paper is to turn the logic around and use 
features to probe the background scale factor. In order to do this, it is important that we 
classify which type of features are observationally sensitive to different scale factors, so can 
be used to distinguish different paradigms; and which are not. This is also one of the issues 
that we will investigate in the next two sections. 

5 Sinusoidal running as trigger 

We now study the correlation functions caused by the sharp feature around to- As em- 
phasized, we concentrate on the universal behavior of the BD vacuum. For the kinematic 
Hamiltonian (14. 5p . the BD vacuum behavior for the mode function 



(CkiCkaCka) = S{ki,k2,k3) 



1 



(4.14) 



(fci/EsA^a) 
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is the same as fl2.3l) except for different normalization factors that vary much slower than 
the vacuum oscillations. Namely, 



ay Aek 

Although the sharp features involve both the horizon and sub-horizon scale physics, to see 
the most important universal feature, it is enough that we look at the subhorizon behavior 

o. 

Due to the sharp feature, e receives some small, but sudden, change. How e evolves 
afterwards is model-dependent. For example, in inflation, it will approach again to an 
attractor solution in a few Hubble time. But to see the universal effect of the sharp feature, let 
us only focus on this sudden change. The properties of the power spectrum and bispectrum 
depend on the relationship between the sharpness of this sudden change, which we denote 
as Ats, and the mode fcj. 

For power spectrum, if 2ki ^ Ar^^^, the feature is very sharp compared to the oscillation 
time-scale in BD vacuum (15. 2p . We can approximate the change in e as a step function, 

Ae ^ esd{T - To) . (5.3) 

Plugging (Q and into fliTTT]) . we get 

^ ~ - l-cos2A;iro . 5.4 

The running behavior cos(2/i;iro) remains similar if 2ki ~ Ar^^^. For 2ki ^ Ar~^, how- 
ever, the oscillation in BD vacuum is much faster and the change in e is averaged out, so 
AP^/P^o 0. The important point is that, unlike the resonance case, the sinusoidal running 
behavior is not unique for inflation, but universal for arbitrary time-dependent background. 
The case for bispectrum is similar. The sinusoidal running for S is universal, 

S ~ /atl cos{Ktq + phase) , K = ki + k2 + k^ , (5.5) 

although the amplitude /tvl now depends on the behavior of the mode function after horizon- 
exit, which is highly model-dependent. If we restrict to inflation, a similar estimate as in 
the power spectrum case can be made for the bispectrum. If ^ Ar~^, we use the 
approximation (15. 3 p and get 

where oq is the scale factor at tq. If K ~ Ar^^^, the sharpness of i) and the oscillation time 
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scale in fl5.2l) are comparable, so fl4.13p leads to 



If ir» Ar.-i, Jnl^O. 

So the correction to the power spectrum is generally very small. For example, specializing 
the model of Sec. IHto the slow-roll inflation case, we have e^/e ~ /3, where /3 is the fraction 
of the kinetic energy of converted to that of a. The size of the bispectrum (15. 7p depends 
on Ats and can be very large if the feature is sharp [HaoArs = HAtg — )■ 0). 

To summarize, first, the most important property caused by sharp feature is the sinusoidal 
running for modes K ^ Tq"^; this shows up as the correction to the power spectrum, 

APr 

— — oc sin(2A;iro + phase) , (5.8) 
no 

and as the leading contribution in the bispectrum, 

5* (X sm.{KTQ + phase) . (5.9) 

Second, the starting point for this running is around the scale fco = |to|^^, which is the 
mode that is crossing the event-horizon at the time of the feature to! the wavelength of 
this sinusoidal running in 2A;i-(or i^'-)space is given by the same scale 27r/co. Third, this 
qualitative behavior is universal for arbitrary time-dependent background (i.e. for all values 
ofp)9 they cannot be used to distinguish inflation from the alternatives, but can be used 
to identify the location of the sharp feature, which is a signal that some massive fields are 
likely to be excited. In the next section, we will study the effects of these massive fields on 
density perturbations, including their profiles, locations and magnitudes. These will become 
the main signals we use to distinguish different primordial universe paradigms. 

We briefly comment that there may be other model-dependent signatures due to inter- 
action with massive modes. For example, if we consider the two-field model in Sec. H] in 
inflationary spacetime, during the sharp turn, quantum fluctuations of the massive modes 
are projected to the curvature mode. By matching the curvaton mode function before and 
after tq, we can see that the factional correction to the power spectrum due to this effect is 
APc_/Pc_Q ~ 6'o(m/if)(-A;iro)^/^ cos[(m/if) In/Ei + phase(/ci)], where ki < r^^ and is the 
turning angle. The phase(fci) is a fci-dependent random phase from the massive modes, due 
to which there is no definite prediction on the oscillatory running. But the main signature is 
the overall amplitude with a blue tilt, ~ , because massive fluctuations decay in expand- 



''For inflation, this type of running has been shown for power spectra |32I34||40] and bispectra |30l29ll35ll36] . 

For cases where the feature is very sharp so modes well within the horizon can be affected, the formulae 



()5.7p gives a better estimation for the maximum bispectrum amplitude than those in [30,5 . For example, 
for a small step in slow-roll potential with width d and relative height c, ^ c{c + e)/{(Pe). 
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ing spacetime. These predictions are more model-dependent and we will not discuss them in 
more details in this paper. They may provide supportive evidence on the detailed process. 



6 Resonant running as evidence 



We now compute the resonance effect on power spectrum and bispectrum induced by the 
excited oscillatory massive fields. As in the previous section, we first compute the correlation 
functions in the general power-low background, and point out the most significant general 
behavior. When we wish to see whether such effects are large enough to be observable, 
we use inflation as the explicit example. This is because our main purpose here is to find 
distinctive signatures for inflation. Otherwise, one can use concrete alternative models as 
explicit examples. 

All the necessary ingredients for this computation are ready. The same formalism in 
Sec. m applies here. For power spectrum, we use (13. 5p for Ae in (14. lip . From Sec. |2l we 
know that for resonance we only need the universal BD behavior (15. 2p for Uk in (14. lip . After 
performing the same type of integral encountered in Sec. [2], we get 



Pr 



CO 



4 eM2 



5/2 



2ki 



-3+ 



2p 



sm 



p2 2m^ /2A;i 



i/p 



2a + 



3tt 



(6.1) 



where Hq are evaluated at to, and, for power-law scale factor, e = 1/p is constant. For 
bispectrum, we use (13. 6p for rj in (I4.13p . The general amplitude is mo del- dependent, but the 
resonant running behavior is given by 



S oc 



-3+f 

2p 



sm 



2mo- 
-p Ho 



i/p 



phase 



(6.2) 



where we have also included a i^'-dependent modulation factor which typically arises but is 
not as robust as the rest of the running behavior. 

In these results, we have defined a parameter kr which denotes the first i^'-mode that 
resonates as soon as the massive field starts to oscillate at to. Namely, kr = 2m„ao- Recall 
that, at to; the comoving mass-scale of the event horizon is ko = |to|~^; and ko is the starting 
mode in the sinusoidal running due to sharp feature. It is important to notice that there is 
a relation between the ratio kr/ko and the ratio 2mfj/Ho, 



\p\ 2m„ 
\l-p\ Ho 



(6.3) 



We can also qualitatively understand how the resonant running is capable of recording 
the scale factor evolution. The oscillating massive field provides periodically oscillating 
background, as well as a resonance scale with constant physical wave-number. Different K- 
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modes of the BD vacuum are stretched or contracted by the scale factor a(t), and resonate 
with the background when their physical wave-number coincide with the resonance scale. 
When the change in K corresponds to the change in t that is equal to the oscillation period of 
the massive mode, the final phase grows by 2n. This is why the arguments of the sinusoidal 
functions in ( ]6.ip and ( 16. 2p are power-law function with the inverse power 1/p, which is the 
inverse function of the power-law in scale factor. 

Take the exponential inflation limit, p ^ 1, in the two-field modelj^ For power spectrum, 
we get 



P. 



CO 



4 eM| V H 



-3 



sm 



2m 



H 



- In 2ki + a 



(6.4) 



where the phase a = {2m„/ H){1 — ln2mo-) + 2a + 7r/4. For bispectrum. 



S 



/m,\9/2 fK 



eMl V H 



sm 



2m^ 
H 



Ini^ + a 



(6.5) 



where a = (2m^//7)(l -ln2m^) + 2a-37r/4. Both (El]) and (ES]) take the CEL form ( 12^]) . 



6.1 Resonant running 

The resonant running for different p are very different. Let us look at the details in the 
power spectra. The bispectra are same after replacing 2ki with K. 

First, kr is the first resonant mode, but for different p the subsequent resonant modes 
are different. For the expanding background p > 1 and p < 0, both (16. ip and (16. 2p apply 
only for 2ki > k^, since lower fc-modes resonate earlier. For the contracting background 
< p < 1, the situation is opposite and the results apply only for 2ki < k^. 

Second, if we denote the local periodicity of the resonant running in fc-space as A/ci, 
we have A/ci oc k^ Qq for p > 1 and p < 0, Afci increases as ki increases; while for 

< p < 1, A/ci increases as ki decreases. 

Several examples of resonant running are plotted in Fig. [31 As emphasized, the dif- 
ferences in these resonant runnings are kept intact even after general curvaton-isocurvaton 
transformation, and they are the faithful signals we can use to distinguish the primordial 
universe paradigms. 



6.2 Running of amplitudes 

Besides the difference in the resonance running, different spacetime backgrounds also give 
rise to different scale-dependence in the modulation amplitudes. These scale dependence 

^For inflation, the effect of the resonance mechanism on power spectrum and non-Gaussianity due to 
periodic features in single field models is studied in t29,i37,.31ji38^ ; the effect of oscillating massive field at 
the beginning of infiation on power spectrum is studied in |39| by introducing a direct coupling to inflaton; 
the effect on power spectrum after integrating out the massive modes is studied in j41| . 
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Figure 3: Resonance running in different time-dependent backgrounds due to features peri- 
odic in time, sin 



Y::^C{K/krY^^ + phase . Note that this does not include the running of 

the amphtudes. In these plots we use C = 21710-/ Hq = 50, kr = 100, phase = 7r/4; and from 
top to bottom, p = 10 (inflation), 2/3 (matter contraction), 0.3 (Ekpyrosis). 
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are much milder because the range of scales Aif over which the variation takes place is 
larger than the local scale K itself. For more complicated multi-field models, such scale 
dependence can be changed due to scale-dependence in the curvaton-isocurvaton couplings. 
So in general they are not always the faithful signatures for our purpose. However, there are 
some dramatic properties which may be useful as supportive evidences, so let us nonetheless 
examine these properties. We define the running indices 



for power spectrum, and 



d\nk ^ ' 



for bispectrum. The (AP^/P^o)a and /nl denote the modulation amplitudes, i.e. the overall 
factors in front of the sinusoidal functions, in AP^/P^o and S, respectivelyjj So for power 
spectrum we have 

rip = -3 + |- , (6.8) 
zp 

and for bispectrum we typically have 

nb=-S + ^ . (6.9) 

For exponential inflation, p ^ 1, both indices are red. Hp = Ub = —3. The factor —3 is 
present for all expanding backgrounds because the amplitude of the massive mode is damped 
by the expansion as t"^^/^. This factor is also present for the contracting backgrounds for 
the following reason. For contracting backgrounds (0 < p < 1), the amplitude of massive 
mode is growing as runs from — oo to in this case). But an important 

difference between the expanding and contracting background is that, in the former, smaller 
/c-modes resonate earlier, but in the latter, larger /c-modes resonate earlier. This is why 
although the amphtude of the massive mode evolves oppositely in time for the two cases, 
their contribution to the running index turns out to be the same. 

There is also an additional factor ~ For contracting background with small p, this 
makes both indices blue, Up ^ 5/{2p) and rih ~ 7/{2p). This factor is due to two reasons. 
First, for contracting background, the resonance scale is fixed while the event horizon is 
shrinking. So the resonance strength (oc |t|^^^) gets weaker for smaller fc-modes. In the 
meanwhile, the couplings (13. 5 p and (13. 6 p depend on H, which (oc |tp and \t\^ respectively) 



^Note that for bispectra we have separated the resonant running from the definition of /nli a-s in [5 . So 
the index is slightly different from the definition unq — 1 in j421l43j . where it is defined for bispectra with 
non-oscillatory running. 
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also get weaker for smaller k. 

As discussed, the details of the indices may not be faithfully kept in terms of the curvature 
mode in more general models; but the dramatic difference between the different cases, such 
as p 3> 1 and < p ^ 1, can serve as supportive evidence. 

We also comment that, because the blue or red running indices are generally of order one 
or larger, the signals we are looking for typically decay away in a few efolds. But this does 
not limit their usage. As we can see from the last subsection and Fig. 121 the differences in 
resonant running for different p are already very clear within a couple of efolds. 



6.3 Amplitudes 

We use the two-field model (14. ip and f l4.2p to show that the amplitudes of the resonant power 
spectra and bispectra can be easily made very large, at least for the inflation models. 

Consider the example of slow-roll inflation. We denote the fraction of the kinetic energy of 
0, that is converted to the energy in the cx-field during the turning and induces its oscillation, 
as 13. So mla\ ~ l3(j)'^. Also note e ~ ^"^/{M^H^). From ([S3D and ([S3]), we have 

fe), - Kir) • 

hL. ~ p(^y" . (6.11) 

So even for a tiny fraction of energy transfer, the resonance amplitudes can be quite large. 
For example, for (3 ~ 10~^, m„/H ~ 10^, we have AP(^/P(^o ~ 0.1 and f^i ~ 10^. Because 
of the spatial inhomogeneity characterized by 5t ~ 10~^/if, the zero-mode oscillation in the 
classical background receives a random phase correction, uj5t. This phase has to be much 
smaller than 27r so that the signals we are interested are not averaged away. Therefore we 
can at most explore the massive modes over five order of magnitudes above H, ma/ H < 10^. 



7 The signature pattern for inflation 

We provide a summary on the signature pattern for the inflation paradigm. Similar summary 
can also be done for each alternative paradigm by specializing the previous general results. 

We have shown that a detection of the resonant form of CEL type induced by massive field 
oscillation in power spectrum or non-Gaussianities is an evidence for the inflation paradigm. 
But to be unambiguous, it is important to strengthen the evidence that this is due to the 
periodically oscillating massive fields, by using other characteristic properties besides the 
resonant running. The following are the signature pattern that we can look for in density 
perturbations: 

• Trigger. Observational signatures associated with sharp feature can be used as a sign 
that some massive fields may be excited. These signatures appear as sinusoidal running 
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in density perturbations, as corrections to power spectrum or dominant components in 
non-Gaussianities. These oscillations start at a scale ko, have constant wave-length ~ 
27rfco in 2fci-(or i^'-)space, and propagate towards larger /c-modes with model-dependent 
growing and then decaying amplitudes. 

• Signal. Excited massive field induces highly oscillatory resonant running in density 
perturbations, as corrections to power spectrum or dominant components in non- 
Gaussianities. This resonant running has a distinct GEL form (12. 7p . For example 
for bispectra, it starts at a scale kr and propagate towards larger iT-modes, typically 
with decaying amplitude and lasting for no more than a few efolds. The oscillating 
wavelength AK is always smaller than the local K, with a fixed ratio that is deter- 
mined by the parameter 2m„/H, i.e. AK/K = TrH/m^. The starting place k^ for the 
previous sinusoidal running and k^ for this resonant running is related by the same 
parameter 2m„/H., i.e. k^/k^ = 2mo-/if. 

It is also likely that several massive fields with different mass are excited at the same 
time. So we may look for different GEL forms in modes much larger ko, each satisfying 
the relation hn/ko = 2mcri/H. The relation between the different GEL forms can 
also be used to conclude that they are induced by the same sharp feature, even in case 
where the observational signatures from the sharp feature is too weak to be observable. 
Namely, by measuring kr and 2m(j/H for each form, they should satisfy 

^rl _ kr2 _ (7 ^\ 

2m„i/H ~ 2m,2/H ~"' ' ^ ' ^ 

• Caveats and solutions. It is also important to note several caveats and possible solu- 
tions. 

In the inflation case we considered above, the mass may be time-dependent. But such 
dependence has to be very dramatic [m^/{rri^H) > 0{1)] to make the flnal resonance 
form differ signiflcantly from the GEL formO So for inflation with massive modes, the 
GEL form is the generic form we expect to measure. The question we concern is how 
non-inflationary spacetime may produce the same speciflc pattern. 

For non-inflationary spacetime, periodic oscillations from massive modes generate dif- 
ferent types of resonant forms (12. 6p . So to mimic the GEL form we need to engineer 
artiflcial features. As we have shown in Sec. [21 features that introduce a background 
oscillation component of the form B{t) ~ e'9Ht/to) can also induce the GEL form. To 
reproduce the signature pattern for inflation, we need to place a sharp feature right at 
the beginning of these repeated features, to satisfy the relation for kr/ko. A possible so- 
lution to such an ambiguity is to detect or constrain more observables, which naturally 

^The dramatic time-dependence in mass will also lead to large running in the oscillating amplitudes (|6.10|) 
and (|6.11|) , therefore modifying the overall running behavior of the amplitudes. 
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arise in the same inflation model, but makes reverse engineering in the alternatives 
more artificial. For example, as we mentioned, it is natural that more than one set of 
CEL forms are present with different 2ma/H. To engineer them in non-inflationary 
spacetime, we need to superimpose repeated features with different g parameter on top 
of each other, and right after the sharp feature. In addition, the characteristic running 
amplitudes of the resonance forms, (16. 8 p and (16. 9p . can be used as supportive evidence. 
Furthermore, a sharp feature in non-inflationary case is likely to excite massive flelds, 
which induce different types of resonant forms. Constraining these forms can provide 
additional supportive evidence. 

Finally, although we expect small excitations of some massive flelds exist generically, 
they are not always observable. For example, for density perturbations at £ ~ (9(10^), 
the highest mass we can possibly detect through this method is 0{1Q^)H. This is most 
likely to be further limited by experimental sensitivities and sky coverage. Therefore 
experiments that target on high multipoles are most useful for our purpose. 

Overall, like the tensor modes, resonance phenomenon induced by massive flelds has its 
generic and distinctive set of predictions for general inflation models. In addition, the 
signatures for different paradigms are different and can be used to distinguish inflation 
from the other paradigms without degeneracy; this aspect is even more advantageous 
than the tensor modes. But it also has similar caveats. Not all parameter space 
are measurable. Also they may be engineered in alternative paradigms by different 
processes, but such engineering can become highly artiflcial by predicting, constraining 
and measuring more observables naturally present in such phenomena. 

8 Experiments and data analyses 

As we know, the tensor mode is determined by the horizon mass-scale in the primordial 
universe, such as the Hubble parameter H in inflation, which may be much higher than 
energy scales accessible in accelerators. Interestingly the mechanisms studied here involve 
energies much larger than H , and therefore is a probe of even higher energy scales. In this 
paper we have shown that such mechanisms can record the time-dependence of the scale 
factor of the primordial universe in terms of distinctive oscillatory running of resonance 
forms. Such features are determined by the properties of the BD vacuum that is shared by 
all scenarios, and are kept intact for general multifleld models. 

In this section, we discuss several experimental and data analyses aspects. As indicated 
by the CEL form, such effects show up in terms of oscillatory signals in /c-space. To observe 
them, the binning in the multipole space A£ has to be much smaller than ^ itself. This 
requires high precision experiments capable of observing density perturbations for large ^. 
The Planck satellite is observing the CMB at maximum multipoles of a few thousands. 
The ground-based telescopes, Atacama Cosmology Telescope (ACT) [H] and South Pole 
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Telescope (SPT) [U], can go up to ten thousands. More speculatively, the 21cm hydrogen 
line may be observed in much lower redshift and in much higher multipoles. 

The CEL form, and most other resonance forms, have highly oscillatory and characteris- 
tic running behavior. This makes it very difficult for other effects, such as the astrophysical, 
nonlinear gravity and systematic effects, to mimic such signals. For example for CMB, the 
observational sensitivities for conventional power spectrum and bispectra are dramatically 
reduced as we go to high £ of several thousands, due to astrophysical effects such as the point 
sources and the Sunyaev-Zeldovich (SZ) effect. This is the case for the ACT and SPT ex- 
periments in the range ^ > 2500. Since the CEL form is orthogonal to these contaminations, 
part of this range may now become important in terms of probing the primordial cosmology. 
For CMB experiments, the nonlinear effects in CMB evolution, which limit the sensitivity 
for various scale-invariant bispectra to be of order /^r^, ~ ^^(l), are also orthogonal to the 
type of signals we study here. Therefore a new assessment is necessary to find out the main 
limiting factors and make forecasts. 

To search for such signals in CMB, instead of starting with a specific template, we need to 
scan a variety of non-separable functional forms. The modal decomposition method [l^HT] 
developed by Fergusson, Shellard and collaborators seems ideal for such goals. This method 
has been mainly applied to general bispectra |^H9] and trispectra [50|[5T| with less dramatic 
scale dependence, but should be able to be generalized to cases with highly oscillatory scale 
dependence, as well as to the power spectrum. 
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